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INTERPRETATION OF COMET SPECTRA 
C. Arpigny 

INTRODUCTION 

This report will be devoted to a discussion of the theoretical 
interpretation of the spectra of comets and in this discussion we 
shall consider successively a number of molecules that have been 
studied recently : CN , CH, C2, C^, OH, CH^ . The first two of this 
list, CN and CH, have been analjsed in greatest detail until now, 
so that we shall be rather brief on the other molecules. We shall 
try to indicate what kind of information can be derived from these 
studies concerning the conditions prevailing in cometary atmospheres 
and to show which fundamental data are needed in the calculation of 
the theoretical spectra. 

Before going into the main part of the report, however, I 
should like to present a resume of a paper of a general character 
which is concerned with the Spectral Classification of Comets and 
which is due to Dr. Bouska. 

Up to now a few hundred spectra of several tens of comets have 
been obtained at different observatories over the world. On the one 
hand, there are high-dispersion spectra photographed with slit 
spectrographs attached to large telescopes, and on the other, small- 
dispersion spectra obtained by means of objective prisms. 

It seems therefore useful to introduce a classification of the 
spectra of the cometary heads. It is evident that it is impossible 
to put into practice a classification similar to that used for 
stellar spectra. The classification of cometary spectra could rather 
be like the classification of meteoric spectra. 
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The cometary spectra show mostly two components : 

(C) continuum which is the solar spectrum reflected on the dust 
particles present in the cometary coma, and 
(E) emission bands (CN, C 2 > and some others) connected with the 
intrinsic radiation of molecules in the cometary atmosphere. 
The apparent intensity of these components may be weak (1), 
normal (2) or strong (3). In some exceptional cases continuum or 
emission bands may be absent ( 0 ), Consequently the classification 
may be as follows : 


Continuum 

Emission Bands 

Apparent Intensity 

CO 

EO 

absent 

Cl 

El 

weak 

C2 

E2 

normal 

C3 

E3 

s trong 


In most cometary spectra the CN (0,0) band is dominant, in 
others the Swan bands of C 2 ; this may be expressed by the following 
symb o 1 s 

c - cyanogen bands dominant, 
s - Swan bands dominant. 

The presence of the sodium doublet D observed in cometary 

I » z 

spectra at heliocentric distances r ^ 1 a.u.may be denoted by the 
symbol n. 

Metallic lines are observed exceptionally, for instance in 
the spectra of Sun-grazing comets. Such spectra may be denoted by 
the symbol M. 

The cometary spectrum depends, of course, on the comet's helio 
centric distance and this distance must therefore be given in the 
classification of the cometary spectra. 
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The following examples illustrate the proposed classification : 

C3E2c(0.6) Comet Ar end-Ro 1 and (1957 III). Strong continuum, normal 

intensity of emission bands, CN (0.0) dominating; r = 0.6 a.u. 
ClE2s(1.2) Comet Ikey a-Everhar t (1966 IV). Weak continuum, normal 

intensity of emission bands. Swan bands C 2 (Av = +1) dominating; 
r = 1.2 a.u. 

I should like to make two very quick comments to this proposed 
classification scheme, which would indeed be useful for rapid 
communication through the lAU circulars, as well as for statistical 
studies of comet spectra. 

First, the strength of the continuum, especially with respect 
to the discrete emissions, depends strongly upon the dispersion of 
the spectrogram used, so that it would seem desirable to specify 
this dispersion in the classification symbol. This point is 
illustrated on Figures 1 and 2, which show the region of the CN 
violet (0,0) band in Comet Bennett (1970 II) as observed at r = 1.0 

O O 

a.u. at the Haute- Provence Observatory using 12 A/mm and 7 A/mm 
reciprocal dispersions respectively. The continuum, which is strong 
even in this ultraviolet part of the spectrum, is nevertheless 
appreciably weaker relative to the discrete emissions when recorded 
at a higher dispersion. 

The second remark is that the classification would be even more 
useful if it contained some information about the spectra of the 
ions, in particular of CO^ . 

Ishould now apologize to those who are well acquainted with 
cometary spectroscopy, but it seems worthwhile, for those who are 
not so familiar with this subject, to start with a short summary of 
the basic ideas underlying the interpretation of the spectra of 
comets. This anyhow will set the stage for what we are going to 
discuss later on. A more complete presentation can be found elsewhere 
(Arpigny 1971). 

When one examines all the possible mechanisms he can think of 
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Fig. 1 - The region of the CN Violet (0,0) band in Comet Bennett 

(1970II) - (Haute- Provence Observatory; original reciprocal 

O 

dispersion : 1.2 A /mm) . 



Fig. 2 -I- The region of the CN Violet (0,0) band in Comet Bennett 

(1970II) (Haute -Provence Observatory; original reciprocal 

O 

dispersion : 7 A/mm) . 
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for exciting the electronic transitions observed in cometary spectra, 

he finds that the only plausible process is one of resonance- 

fluorescence in which solar radiation is absorbed in their resonance 

transitions by the cometary molecules, which re-emit right away in 

the same electronic systems . This is due to the very long timescales 

associated with the other processes, typically of the order of 10^ to 
6 

10 sec or larger. The processes in question are, for instance, 
collisional excitation by cometary heavy particles or electrons or 
by solar wind particles and the long characteristic times are due 
either to low densities, or to the low energies available, to rather 
small cross-sections, or to a combination of these. By contrast, the 
characteristic times for absorption of solar radiation in the "optical" 
region are of the order of 10 to 100 sec for the molecules or atoms 
observed in comets near one a.u. from the sun. The corresponding 
relaxation time or the time for a steady state to be established is 
about several 100 to several 1000 sec. Now the excitation conditions, 
which depend upon the intensity of the solar radiation available at 
the various transitions and hence upon the heliocentric distance r 
and upon the heliocentric radial velocity dr/dt, remain essentially 
constant during a time of that order and we are therefore justified 
in assuming that stationary conditions will indeed be therefore 
justified in assuming that stationary conditions will indeed be 
reached . 

CN. 

As an example, we shall consider the emissions of cyanogen and 
we shall first simplify the situation in order to emphasize the 
essential features of the mechanism. Thus, we shall (a) assume that 
the only transitions to be considered are those in the (0,0) Violet 
band (b) neglect the small spin splitting in the X and B E 
electronic states (c) consider that all the CN radicals are moving 
at the same velocity relative to the sun, i.e. the orbital dr/dt. 
According to our assumptions, we have only two kinds of electronic 
transitions, P and R transitions. In addition to these, pure 
rotational transitions in the ground state must be included. Under 
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s t eady- s t a te conditions we equate the rates of population and of 
depopulation of the rotational levels and by considering the 
absorption-emission cycles which connect these levels it is easy 
to obtain at once a set of equations involving the populations, Xj^ , 
of the lower rotational levels only. For level N ^N" , rotational 
quantum number) we have 

^N-2‘ ^00 ^N+ 1 ^+I,N ^N+ 2 ' ^ \+ 2^ ® N^OO 

[ (Sr)n(Sp>j^^2 * ^®R^N-2^ ^00 

where the C's and the A's are transition rates per molecule in the 
initial level; Cqq is the overall absorption rate in the (0,0) band, 
s and s are the absorption and emission intensity factors, respectively, 
representing the appropriate rotational line strength divided by the 
statistical weight of the initial level of the transition, s contains 
in addition the residual intensity factor, i^^, which takes account 
of the energy available in the solar spectrum at the excitation 
wavelength of the transition. For example, the first term in 
equation (1) is the rate of the RP fluorescence process N-2 N ' = 

N-1 N (x ^_2 (Sp^) Cqq = transition rate for absorption of the R 
(N-2) line; “ branching ratio for the downward P(N) transition). 

Once the system represented by eq. (1) has been solved (in 
practice about 30 rotational levels must be included in the lower 
electronic state), the populations of the upper levels as well as 
the relative intensities of the lines can be readily computed. These 
populations and intensities will be governed by two effects : 

(i) the first effect is the balance between the two kinds of processei: 
a) the fluorescence cycles tend to increase the rotational 

energy of the molecules in the electronic ground state (if 
we start out with all molecules in the lowest rotatixial level, 
the distribution of molecules will be spread toward higher 
levels via RP cycles) ; 
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b) the pure rotational transitions on the contrary tend to 
return molecules back to lower levels; 

(ii) besides this competition, the spectral energy distribution in 
the exciting solar radiation plays a predominant role. 

The fluorescence terms in the steady-state equations all contain 
while the pure rotation transition terms are proportional to 
Arot and the determining factor is the ratio between these, represented 
by a quantity 


f 

oo 

y being the permanent dipole moment of the molecule, f^^ the f-value 
of the electronic transition. The distribution of populations will 
go through a maximum which will occur at a smaller rotational energy, 
the larger (R is. Thus, for a given molecule, with y and f^^ fixed, 
the maximum will occur at lower rotational quantum numbers for the 
larger heliocentric distances. On the other hand, at a given r, the 
rotational "excitation temperature" will be different for the 
different molecules. It will be low for heteronuclear molecules, and 
large for homonuclear molecules like C 2 or for an atom like which 
has a number of low-lying metastable terms. In the latter case ( (R 
very small) it can be shown that the excitation temperature 
approximates the colour temperature of the exciting solar radiation 
in the relevant wavelength range. 

If we now consider the second effect mentioned above, i.e. the 
effect of the solar spectrum, we see that the very irregular 
character of this spectrum, with numerous Fraunhofer lines in the 
region of the CN Violet system for example, is responsible for the 
irregularities observed in the intensity profiles of the cometary 
emissions. This interpretation was given for the first time some 
30 years ago by Swings (1941) who also pointed out that the intensity 
distribution would be very sensitive to the heliocentric radial 
velocity of the comet. 
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The first computations based on such a simplified ( "zeroth order 
approximation") model were made about 10 years ago by Carrington 
(1962) and by myself (1964). The comparison of the observed and 
theoretical spectra for Comet Mrkos (1957V) showed that the general 
character of the distribution was rather well represented by the 
computations, but there were some discrepancies. In particuliar, 
the R(7) and R(8) lines were predicted appreciably too high. It 
was hoped that the situation would be improved if one included 
transitions in the Red system of CN in the calculations because 
there are fewer absorption lines in the red region of the solar 
spectrum and this would smooth out the distribution of populations. 
While this is indeed the case, one is still left with differences 
which can reach about 35 % for a few weak lines in the region of 
higher quantum numbers. Another example of simplified model 
calculations is given in Figure 3. This is the R branch of the CN 
(0,0) band in Comet Bennett (1970 II). The spectrum was computed 
by Aikman et al (1974) and one sees again a good general agreement, 
although a few significant differences are noticed between observed 
and theoretical spectra. 

Let us note that from the value of (R required to compute the 
synthetic spectra one can derive a value for the dipole moment if 
f^^ is known and we may mention in passing that the value obtained 
from the zeroth order approximation was used in the interpretation 
of the excitation of interstellar CN radicals (Field and Hitchcock 
1966) at a time when no experimental determination existed for y. 

This is an example where, even if the value for y was crude, cometary 
spectroscopy has been of some help in connection with an outstanding 
as trop hy s i cal problem, which in fact was related with no less than 
the origin of the Universe ! 

We shall come back to the analysis of the rotational structure 
of CN, but let us first say a few words about the relative band 

intensities in the Red and Violet systems of this radical in comets. 

2 + 2 

We consider nine vibratmal levels in each of the X Z , A II, and 
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The R-branch of the CN Violet (0,0) band in Comet Bennett 
(1970II). Dashed lines : observed spectrum; solid lines 
spectrum calculated on the basis of the simplified model 
assumptions (taken from Aikman et al , 1974). 





B 2E+ electronic states and write steady-state equations under 

solar radiation excitation. Just as we have defined a quantity (R 

for the rotational structure, we define here v = a'^*^/c'^ 

10 00/ 10 

being the Einstein coefficient for spontaneous emission of the 
vibration-rotation fundamental band and the absorption rate in 

V i b 

the Violet (0,0) band. Ajq is not known, so that v must be treated 

as a variable. Another relevant quantity is the ratio of the f- 

R 

values of the Violet and Red (0,0) bands, (hereafter called 

FR) , about which there is also some uncertainty. Some computations 
have been made recently by T. Banks and myself (1973) in an attempt 
to interpret O’Dell's and some older observations by Dufay and 
Swings of some CN Red and Violet bands in comets Tago -Sa t o -Kos ak a 
(1969 XI) and Mr kos (1957 V), respectively. By varying v (or v its 
value at r = 1 a.u.) and FR and looking for a fit with the observed 
intensity ratios, we find a region defined in Figure 4 by the thick- 
lined box inside which the fit is possible. We have here two cases 
corresponding to two different sets of relative f-values within the 
Red and Violet systems. The dot and error bar indicate the most 
recent experimental value for FR and we see that in one case there 
is no overlap while in the other thffe is. The thin lines indicate 
how the right-hand boundary of the "fitting" box is shifted for a 
change by 5 % and by 10 % in the observed value of ^3j/^20’ 
ratio of the Red (3, 1) and (2, 0) bands. The sensitivity is very 

highindeed. The conclusion of our study is that v is certainly 

^ V i b ” 1 ^ 

> 10, i.e. Ajq > 1 sec , and that FR ^ 15, but new, accurate 

measurements of several band intensity ratios in comets would be 

most welcome. 

We now return to the rotational structure of the CN comet 
spectra. A number of important refinements were brought about by 
Malaise (1970) a few years ago, especially with respect to the 
model of the coma. As far as the CN molecule itself is concerned. 
Malaise considers radiative transitions in the (0, 0) Violet band 
and in several Red bands, as well as pure rotation and vibration- 
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rotation transitions in the ground state. He also takes account of 
the spin splitting, as required. Besides these radiative transitions. 
Malaise introduces for the first time in such calculations the 
effects of some collisional processes which could excite the 
rotational levels within the ground electronic state. Such transitions 
among rotational levels require much less energy and have larger 
cross-sections than the electronic excitations and if the densities 
of exciting molecules are high enough, the collisions may compete 
with or even dominate the radiative transitions within the ground 
state. One immediately realizes the considerable interest of this, 
because it gives a hope to estimate total particle densities or the 
densities of the main constituents of the cometary gas, as opposed 
to the already known densities of the observed radicals, which are 
minor constituents. If the collisions are predominant, the lower 
rotational levels of CN will be populated according to Boltzmann's 
law at the temperature of the colliding particles and the way in 
which Malaise takes account of the collisions is as follows. The 
population of each lower rotational level is written as a linear 
combination of the corresponding Boltzmann population and pure 
fluorescence population : 

X = ax^ + (1 - a)x^, {?.) 


where a gives a measure of the relative importance of the collisional 
and radiative mechanisms and is taken from the relation 


g 

1 -g 



n V a . 


(3) 


T being the time necessary to reach a steady-state under fluorescence 
F 

conditions (this is estimated to be a few 100 sec for CN at 1 a.u. 
from the sun), the collision time, n the density of exciting 

molecules, v the average relative velocity, and o the cross-section 
for collisional excitation. Now this looks reasonable at first sight, 
since the expression (2) gives correct asymptotic values at the 
extremes of high and low densities. When n is high, the collisions 
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dominate , the ratio (3) will be large, hence a will be near unity 
and the populations will be Bo 1 1 zmannian , while for low n the ratio 
is small, hence a will be near zero and we have essentially pure 
fluorescence populations. However, the validity of eq. (2) in the 
intermediate situations is highly questionable, as we shall see 
later on when we discuss the case the of the CH radical. At any 
rate, n and therefore a are functions of position within the comet 
and the calculations have to be made and eq. (2) used at a number 
of points. An integration is then performed along the line of sight. 
For this, of course, model distributions must be appropriately 
chosen for n (R) and n^^^ (R)» R being the distance from the centre 
of the comet. The former distribution is adjusted while making the 
computations of synthetic spectra, whereas the latter is derived 
from the observed radial profile of CN emissions (intensity 
distribution in the direction perpendicular to the dispersion). 
Furthermore, it can be shown that the cometary gas is optically thin 
in the light of CN, so that the integration reduces to a simple 
summation. On the other hand, the kinetic temperature of the exciting 
molecules (related to v) is considered to be constant throughout the 
comet, which is only a rough approximation. 

Another important effect comes in also, that is the influence 
of the motions of the radicals within the atmosphere of the comet. 

Such motions are of course expected, but their existence is actually 
inferred in the interpretation of the so-called Greenstein effect, 
which consists in differences in the relative intensities of the 

J 

spectral lines in different regions of the comet, for instance on 
the two sides of the nucleus. The model adopted is one of uniform, 
isotropic expansion where the radicals move radially away from the 
center with velocity v^ . What counts, of course, is the component 
of this v^ on the direction to the sun and Malaise computes a 
projection factor in terms of the phase angle sun-comet-earth, 
assuming that the sun lies in the plane defined by the slit of the 
spectrograph and the line of sight. This assumption is often not 

c.y 


809 



valid, however, and this does introduce in some cases significant 
errors in the residual intensities i which enter the steady-state 
equations. The heliocentric radial component of v^ obviously varies 
with position along the line of sight and this is another reason 
for making an integration. 

Four high-dispersion and four medium-dispersion spectra of 

comets obtained at the Haute-Provence Observatory were studied by 

Malaise. The comparison between observed and theoretical spectra 

shows a rather satisfactory agreement, but there are again a few 

lines for which the discrepancy reaches about 20 %. Two quantities 

come out of this kind of analysis. One is T, the kinetic temperature 

of the colliding molecules and the other is the expression n^OTj,//M, 

where n^ is the total density (or a minimum value for this total 

density, if the exciting species is not a predominant constituent) 

4 

at some reference distance from the nucleus, say at R = 10 km, while 
M is the reduced. mass of CN and the exciting molecule. There are 
three comets for which Malaise finds important collisional effects, 
namely Seki-Lines (1962 III), Ikeya (1963 I), and P/Encke (1961 '!) 

4 

(a^, the value of a at R = 10 km, equal to 0.90, 0.90 and 0.55, 
respectively). Assuming the gas to be made primarily of H 2 O or CO for 
instance, and taking a equal to the gas-kinetic cross-section, one 
derives (after multiplication of the second member of eq. (9) of 

Malaise (1970) by a correcting factor of 1 / tt , that of eq. (11) by 

9 “ 3 8 ” 3 

2/n) : n == 10^ cm for the first two comets, and n =10 cm 

o ° 33 

for P/Encke. The corresponding total production rates are Q = 10 

32 - 1 

and 10 molec, sec , respectively. These are lower limits if the 

production of the exciting molecules does not take place entirely in a region 

4 

the dimensions of which are small compared to 10 km. The first figure 
is some three orders of magnitude larger than Q(H 20 ) observed in 
Comet Bennett (1970II), which belongs to the same general class of 
rather bright objects as Seki-Lines and Ikeya, while the second 
figure exceeds by about four orders of magnitude the value of Q(H 20 ) 
derived for Comet Encke itself at a similar heliocentric distance 
during itsperihelion passage of 1970. One may also say that with 
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the above densities and production rates . the total mass loss kt/U of 
each of these comets during one revolution around the sun would be 
close to or even greater than its mass^y^, as shown in Table 1. The 
AW s computed here are based upon the Q's quoted above, parabolic 
orbits, and various power laws for the variation of Q with heliocentric 
distance Q t r with n = n + n,/r , n and n, being in each case 
chosen in such a way that n lies near 2 to 2.5 for r near one a.u. — 
cf Keller and Lillie, 1974). On the other hand, the^Yi^'s are evaluated 
using the radii derived by Miss Roemer (1965) from "nuclear" 

magnitudes (the albedo is taken here equal to O.I, probably on the 

• . . -3 

low side of the actual value) together with a density of one g.cm 

The corresponding values ot represent upper limits to the masses 

of the comets, for the nuclear magnitude tends to overestimate the 

luminosity of the nucleus, hence its radius, because of a non- 

negligible contribution from a halo surrounding the nucleus, and if 

we have to do with separated particles rather than with just a 

monolithic nucleus, the total mass for a given total reflecting area 

will be smaller than if we assume the presence of only one single 

body. Clearly there is a very serious difficulty regarding the 

molecular densities which Malaise's results lead to. There are only 

two factors multiptying n^ in the expression given earlier that could 

possibly influence the conclusions, namely o and x„. 

- 1 5 2 

The value adopted for a, i.e. 2 x 10 cm , does not appear 

to be too small by any significant factor and on the contrary can be 

deemed rather good from what is known about cross-sections for 

translational to rotational energy transfer (Gordon et al, 1968). One 

might consider that the agents producing rotational excitation are 

charged particles. With the very large cross-section calculated by 

Crawford et al (1969) for excitation of CN by slow (® 0.05 ev) 

electrons, in order to reach the collision rates of 0.03 to 0.3 sec 

5 3 

in question here one would need 10 such electrons per cm or more, 
which seems rather high, although we have unfortunately no direct 
evidence about electrons in comets. As far as heavy ions are concerned, 
while it cannot be expected that the corresponding cross-sections be 
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Table 1. 

A Comparison Between Mass Losses and Masses for Three Comets. 


Comet 

^//(lO^^g.) 


Seki-Lines 
( 1 962III) 

>100-300 

<40 

Ikeya 

(19631) 

> 20-30 

<20 

P /Encke 
(19611) 

> 5-7 

< 2 


is calculated as follows. It is assumed, for the sake of 
simplicity in these order of magnitude estimates, that the rate of 
mass loss is the same after perihelion as it was before, so that 

A^/= 2 /°° Q dt (t = o at perihelion), 

o 

Now, if Q is the value of Q derived from a spectrum taken at r = 
rs and if ng = n(rg), we can write 
r 

Q = Q (— )’^s (— )^ (q = perihelion distance), 

S ^ IT 

2 2 

while for a parabolic orbit we have r = q (1+g ) and x = 3+33 » 
provided that g = tan (w/2) (w = true anomaly) and we have 
substituted x = C” ' t , with C = ( 2 9 ^) 1 / 2 . ( 3k) “ 1 (k = Gaussian 
gravitational constant). Finally, we obtain 

= 6CQg(-|)^s /”(1 + S^) 

which is readily evaluated by a numerical integration. Note that 
the use of a parabolic orbit is a good approximation even for 
Comet Encke because the contribution to the integral gets small 
once r becomes large compared to 1 a.u. or 3 large compared with 
(q~l - 1 ) 1/2 (the perihelion distances of the three comets are 

0.03, 0.63, and 0.34 a.u. respectively). Besides, the A^h are 
essentially unchanged if it is assumed that the activity of the 
comets stops at some distance r|v| from the sun, say 3 a.u.^, and 
the upper limit of the integral accordingly replaced by 3 m = 

(rM/ q - 1)1/2. 
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three or four orders of magnitude larger than the gas kinetic cross- 
sections, it is also unlikely that the production rate of ionized 
species be much larger than the production rate for the neutrals. 

On the other hand, one might think that one way out of the 

difficulty would be to say that tp has been wildly underestimated. 

In fact, however, this relaxation time was computed for CN excited 

by sunlight some ten years ago (Arpigny, 1964) and the result was 

2 

of the order of 10 sec at one a.u. from the s'un, which indeed agrees 

with the value adopted by Malaise. Besides, it may be pointed out 

that if one reduced the densities considerably by artificially 

increasing tp, for instance, by several orders of magnitude, then he 

would be faced with another problem, namely that the collision times 

3 

even firther in, say at 10 km from the nucleus, would become large 
compared with the radiative lifetimes of the rotational levels, so 
that the concept of a Boltzmann distribution of populations would become 
meaningless . 

Let us now make a remark illustrating another point, related to 
the fundamental data required for the computation of the synthetic 
spectra. Malaise found in particular for a spectrum of Comet Ikeya 
(19631) that the reason why collisions were needed to reproduce the 
observed spectrum was that fluorescence alone predicted too small 
intensities for the lower rotational lines relative to the higher, 
whereas a better agreement could be obtained by introducing "mixed" 
populations (eq. 2) and adjusting the "rotational temperature" 
properly. However, it can be argued that the relative intensities of 
the lower rotational lines could be increased, under pure fluorescent 
excitation, by increasing the value of(R , i.e. of the dipole moment 
of CN . Indeed, the experimental value of y is about twice the value 
used by Malaise. More should be said about the values of (R and of v 
adopted by this author, but such a discussion would be outside the 
scope of the present review. 

Something should also be said about the wavelengths scales. In 
spite of Malaise's very careful and thoroigh discussion of the data 
on the X's of the CN Violet lines that were available when he made 
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his comp utations, he had to use X's which, we know now from recent 

o 

analyses, were in a number of cases in error by up to 15-20 mA 
(although many lines had better X's). This is significant because 
such displacements in the solar spectrum can produce appreciable 

O 

changes in i,. The new X's are accurate within ~ + 2 mA . Furthermore, 

^ O 

we now have separate wavelengths (AX = 15-30 mA) for the components 
of the spin doublets (down to N = 7), which were unresolved 
previously. On the other hand. Malaise had to make use of solar 
spectra of the center of the disk which were also affected by 

O 

wavelength uncertainties of the order of 5-10 mA. Recordings of the 
solar spectrum, of the integrated light from the disk, with 

O 

considerably more accurate wavelengths (I AX I < 2-3 mA) now exist. Use 
of the new data should hopefully give a satisfactory interpretation 
of the Greenstein effect and thus allow a good determination of the 
velocities of the CN radicals in the coma (The velocities given by 

O 

Malaise correspond to shifts that are smaller than 10 mA, so that 
the wavelengths scales should preferably be known with a rather high 
accuracy ! ) . 

Finally, let us have a look at Figure 5, which shows a very 
beautiful spectrum of Comet Bennett (1970II) taken by Preston with 
the 100-inch telescope on Mount Wilson. At the reciprocal dispersion 

O 

of 4.5 A/mm the P-branch of the CN Violet (0,0) band is completely 
resolved. We also see a number of weak lines in between the principal 
lines : these belong to the (1,1) band and they should be included 

in a refined treatment, especially for example if one wants to make 
a search for isotopic lines, which also fall in between the lines of 
the main band. Notice a somewhat stronger line just shortward of R(14) 
this is a blend of the P(2) line of the (1,1) band with a perturbed 
line. There are a few such perturbed lines in the R and in the P 
branch of the (0,0) band and their analysis in very high-dispersion 
comet spectra could provide some insight on the formation of CN 
radicals in excited states and on the rates of collisional processes. 

In conclusion, it is clear that CN in comets deserves further 
study. The more detailed analysis of the rotational structure of the 
Violet (0,0) band shows a good agreement between computed and observed 
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The CN Violet (0,0) band in Comet Bennett (1970II) (Mount 

O 

Wilson Observatory; original digsersion : 4.5 A/mm). Note 
the weak lines due to the (1,1) band in between the stronger 
features. The emission next to R(14) on the left is partly 
due to a perturbed line. 
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spectra, but we have seen that the zeroth order approximation already 
gave rather good results. This is indeed a rather awkward situation 
and in order to interpret the finer details and thus to assess the 
actual importance of collisional effects, to derive reliable 
information on the densities, temperatures, and velocities within 
cometary atmospheres, if is in fact necessary to treat the excitation 
mechanism rigorously, to use the most precise fundamental data 
available, to improve upon the model of the coma. It would even be 
desirable, in my opinion, to have recourse to photoelectric photometry, 
to avail ourselves of a higher accuracy in the observations themselves, 
just because the margin between the results of the simplified analysis 
and those of the more refined treatment is rather narrow. 


CH . 

We now come to the spectrum of the A 4300 band of the CH radical. 

Figure 6 shows density tracings of this emission in four comets observed 

at the Haute-Provence Observatory (V90, V927, V964) and at the McDonald 

Observatory (8384C). It is seen that the whole band reduces to a few 

rotational lines and that the relative intensities of these lines are 

rather different in different spectra. The comparison between one 

observed spectrum and the corresponding computed spectrum based on a 

zeroth order approximation (Arpigny, 1965) showed a fairly good 

agreement, with differences of 20 to 30 % in the relative intensities 

2 

of one or two lines. However, the value to be adopted for p ^^oo 

that case was about 5 times smaller than the experimental value for 

this ratio. I should now like to present some preliminary results 

obtained very recently by Miss Klutz and myself and based on a 

considerably more elaborate treatment. 

We include fluorescence transition in the (0,0) bands of the 
2 2 2 2 

A A - X n and in the B I - X II systems of CH and take acount of 

the spin-doubling in the various e lect ronic - s t a t e s , as well as of 

2 

the A-doubling in the II state. On the other hand, we adopt the 
isothermal model for the coma, with uniform, isotropic expansion 
with velocity v^. The heliocentric radial component of Vq is computed 
correctly, taking account of the angle between the sun-comet-earth 
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plane and the plane defined by the line of sight and the slit of 
the spectrograph. The colliding particles that could excite the lower 
rotational levels are distributed according to n(R)v R being 

the distance from the centre of the comet and k 2 . For the same 
reasons as before, we have to solve steady-state equat,ions at 
several points along the line of sight and we introduce the effects 
of the collisions right at this stage rather than using linear 
combinations of Boltzmann and pure fluorescence terms for the 
populations of the lower levels. In other words, in addition to 
radiative, we also have collisional matrix elements of the form 
nva . We have , considered several models for the cross-sections 
corresponding to different selection rules and to different laws of 
variation with the rotational quantum number, but we have been unable 
to distinguish between these nodels from this first investigation. 

The results obtained so far will be illustrated by means of 
the example of a spectrum (V927) of Comet Ikeya (19631). We shall 
consider three kinds of situations : 

(1) cases in which the radiative processes dominate in the ground 
electronic state , 

(2) cases where the collisional processes are predominant airiong the 
lower rotational levels , or 

(3) intermediate cases. 

Let us first take the radiation-dominated ( (R -d ) case and refer 

to the upper part of Figure 7, where we plot the computed relative 

intensities of the various lines against the fundamental parameter 

(R j . Then each curve is marked with two vertical bars representing 

the corresponding observed intensity with its error limits. A fit 

with the observations is obtained of course if the various (R j 

intervals so defined do overlap. As we see, this does not happen 

for this (R-d case, whatever (R . may be. The experimental values of 

“63 

u and of f yield a value of 50 x 10 cm for (R,i and it is clear 
that for this value the higher rotational lines are predicted too 
weak, the lower rotational lines too strong. How about the collision- 
dominated ('<g’-d) case ? This is illustrated in the lower part of 
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The relative intensity curves for the principal lines under 
radiation-dominated (top) and collision-dominated conditions 
(bottom). The vertical bars on each curve represent the 
corresponding observed intensity with its error limits. The 
intensities are normalized to R,(l) = 1. 
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Figure 7 where the abscissae are the values of the kinetic temperature 
of the colliding molecules. Here again we find no solution. The value 
of T nearest to a solution would be rather low, near 100°K, but Rj(2) 
and R 2 (3) are mutually exclusive and there is indeed no fit. From 

these extreme cases it can be concnded that n , the density at R = 

4 ... ° 4 -3 

10 km, must be contained in the interval from about 10 cm to 

8 ~ 3 

about 10 cm (this is for gas kinetic cross-sections and H 2 O 
exciting molecules). 

Thus we are led to turn to intermediate cases. Let us take (R j 

equal to the experimental value and fix T at 600°K, for instance. 

In Figure 8 , then, we plot the relative intensities versus log n 

4-3 ° 

and we have a solution, corresponding to n^ = 5 x 1 0 cm . It is 

instructive to examine the spectra themselves, which are represented 

schematically in Figure 9, the intensity of each line being shown by 

a vertical line. The observed spectrum appears as solid lines, the 

theoretical spectrum as dotted lines, all normalized at Rj(l). This 

is done for five values of n^, each value being indicated at the top 

of the corresponding spectrum. As mentioned above, with . = 50 x 
_ 6 — 3 

10 cm the pure fluorescence predicts too weak higher rotational 

lines; there is a deficiency in rotational excitation energy and we 

supply this energy by putting in molecules at some appropriate T. 

3 - 3 

With n = 5 X 10 cm we do not have enough molecules, so that the 

° . . 5-3 

higher lines are still too faint whereas at n^ = 5 x 1 0 cm on the 

contrary these lines become too strong because there are too many 

molecules. It would seem from this that a unique solution exists, 

4-3 

corresponding to n^ == 5 x 1 0 cm . That such is not the case, 

however, is seen when one wonders about the value of T. It turns out 

that situations similar to the one we have just met with occur at 

temperatures in the range of roughly 400 - 800°K. What happens is 

shown in Figure 10. As we go from 600 to 800°K, for example, the 

relative intensity curves given here for two lines move in such a 

way that the n intervals still overlap so that we still get a 

o 

solution. The value of n is well defined in each case because the 
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Relative intensity curves for the principal lines of CH 
under intermediate conditions. The vertical bars on each 
curve represent the corresponding observed intensity with 
its error limits. The intensities are normalized to (R j ( 1 ) 
The arrow indicates the approximate value of log n^ for 
which a fit with the observations is possible. 
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Fig. 9 - Comparison between the observed spectrum (solid lines)and 
corresponding spectra of CH calculated (dotted lines) for 
intermediate cases with various densities of exciting 
molecules. 





overlap of the intervals is small. For T = 400°K,n is found to 

5 -3 . '^2-3 

be = 10 cm , while at T = 800°K it is - 2 x 10^ cm . At lower 

T we need more molecules to supply the missing excitation energy. 

Similar situations are encountered for some other spectra of 

comets Burnham (19601), Ikeya (19631), Bennett (1970II), and 

Kohoutek (1973X11), which we are studying. The , dens i t i e s fall in 
4 6 -3 

the range 10 - 10 cm and the temperatures are ill-defined, 

covering a range of several hundred degrees. Here also it is felt 
that if very accurate observations obtained by photoelectric 
photometry were available, one would be in a better position to 
make the adjustment of the parameters and to improve the model of 
the cometary atmosphere. On the other hand, if the same molecules 
contribute to the excitation of the lower rotational levels of CH 
and CN for example, which seems likely, then the values for the 
densities and temperatures that come out from the analysis of the 
emissions of these (as well as of other) radicals in one given 
comet spectrum should obvious ly agree . This condition will have to 
be borne in mind in future studies. 

We can take advantage of the fact that we are able here to 

calculate steady-state populations with radiative and collisional 

terms and at the same time to compute pure fluorescence populations 

as well as Boltzmann populations, in order to test the validity 

of a linear approximation such as that described in the discussion 

of CN . In other words, referring to eq. (2) we have the means of 

knowing x (corresponding to a given density n of exciting molecules) 
B F 

as well as x and x , hence we can deduce a. If eq. (2) had any 
general meaning, the same value should of course be obtained for a 
when using the populations of all the levels. This has been done 
for a number of values of n corresponding to (R-d, 9f-d, and 
intermediate cases as defined above, with the result that a is 
anything but constant for the various rotational levels. When n is 
very low or very high, a is very nearly equal to zero or to one for 
all levels, as expected. However, it covers a very wide range (e.g. 
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0.10 to 1.45, or 0.20 to 1.25) when n is intermediate. This is not 

unexpected either, because the radiative lifetimes of the different 

lower rotational levels are different (in particular owing to the 

Swings effect), so that the ratio of the characteristic time against 

radiative transition to the characteristic collision time is not 

the same for the different levels. Similar conclusions would most 

probably apply to the case of CN . 

Finally, let us mention that only for one spectrum, namely a 

spectrum of Comet Bennett (1970II) has it been possible to estimate 

the expansion velocity v^ of the radicals, the other spectra being 

rather .insensitive to velocity shifts in the interval 0-1 km. sec 

In the case of this spectrum of Comet Bennett the intensity 

distributions on the two sides of the centre are very different 

from one another and it is found that v is equal to 1.0 + 0.3 km. 

- 1 o - 

sec 



Now a few words about the C 2 radical. The high rotational and 

vibrational excitation temperatures derived from the Swan system 

are interpreted in terms of the resonance-fluorescence mechanism as 

due to the homonuclear character of this molecule, i.e. to the very 

small value of (R However, when one looks at the spectra or at the 

intensity distribution one meets with difficulties here also. For 

example, one finds appreciably different vibrational temperatures 

in different band sequences. Gebel (1970) compares the values 

observed for flux ratios like F(Av = +l)/F(Av = 0) and F(Av = -1)/ 

F(Av = 0) to the corresponding theoretical values for various 

temperatures and finds that the former indicates relatively low 

T ('^ 3000°K), while the latter corresponds to higher vibrational 
vib 

excitation temperatures (> 6000°K) . This can be attributed to an 
underpopulat ion of the upper v" levels with respect to the lower 
v" levels, because the Av = +1 sequence is produced to a considerable 
extent by transitions from the upper v" levels, while the Av =-l 
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sequence is formed mainly from the lower v" levels. No explanation 
has yet been given for this anomaly. An attempt by Gebel to explain 
this anomaly by considering pure vibration transitions, which can 
be shown to be very weak anyway, has failed as can be seen in 
Figure II. I should like to point out in this context that a 
determination of the variation of the transition moment with 
internuclear distance for the Swan system would be most welcome for 
the relative f-values of the various bands are important data that 
come into the interpretation of these C 2 emissions. 

There is alsa the question of the presence of the Phillips 

O 

system. A few bands of this system near 8000 A have been detected 

for the first time by O'Dell (1971) in Comet Tago-Sato-Kosaka (1969IX) and 

it should be interesting to observe them again in various comets 

to derive the singlet/triplet ratio for this might give some 

information concerning the production mechanism of C 2 • 

As far as the rotational lines of the Swan bands are concerned, 

a number of anomalies have been found. For instance, in a detailed 

study of dozens of lines in several high-dispersion spectra 

Woszczyk (1970) has found that seventy percent of the lines could 

be interpreted qualitatively by the resonance-fluorescence process, 

but that the intensities of the remaining lines remained unexplained. 

The study of the C„ emissions has yielded important information, 

^12 13 

namely estimates of the C/ C isotopic ratio in three comets; 

Ikeya (19631) (Stawikowski and Greenstein, 1964), Tago-Sato-Kosaka 
(1969IX) (Owen 1973), and Kohoutek (1973X11) (Banks et al 1974, 

Kikuchi and Okazaki, 1975). The isotopic shift for the shortward- 

O 

degraded (1,0) Swan band amounts to +8 A, so that the isotopic band 
head falls well outside that of the main band. Photoelectric scans 

O 

of this region near 4737 A have been obtained for the first time on 
Comet Kohoutek (1973X11) by Banks, Lambert and myself with the 
McDonald 107-inch reflector and the Tull coude scanner. The scans 
reproduced in Figure 12 have a resolution comparable with that of 

O 

the previous photographic observations (==0.4-0. 5 A). Unfortunately, 
the C C (1,0) band is blended with NH 2 emissions. In the hope 
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Fig, II - Calculated (solid lines) and observed band sequences flux 
ratios for various comets (taken from Gebel 1970). 
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NORMALIZED INTENSITY 



WAVELENGTH 


Fig. 12 - Photoelectric scans of Comet Kohoutek (1973X11) at a 

° 12 13 

resolution of 0.4 A. The C C + NH„ blend is seen at 

O ^ 

4745 A. 
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of separating these emissions we have used a higher resolution 

O 

0.15 A), which is illustrated in Figure 13 where we see that the 

NH„ contribution is made of four lines. Two of these are separated 
■ 12 13 

from the C C band, but the other two are still badly blended with 

this band. Thus the relative intensities of the NH 2 lines must be 

estimated by comparison with the laboratory spectrum and by correcting 

these intensities for the effect of the solar absorption lines in 

an approximate way. In this way we can subtract the contribution of 

NH 2 to the blend. An example of the comparison between observed and 

synthesized spectra is shown in Figure 14. Once we know how much 
12 13 . 

C C contributes to the blend, we derive the isotopic ratio by 
relating the intensity of this blend to that of the main band head. 
Care must be exercised in this analysis because it has been found 
that the overall strength of NH^ is variable. The values obtained 

1213.. ^ 

for the C/ C ratio in the comets quoted above agree within the 
observational errors with the terrestrial ratio of about 90. In 
order to be really significant this isotopic ratio should be 
determined in many comets, using high-resolution spectra as far as 
possible. This unfortunately eliminates a large number of comets, 
which are too faint. Nevertheless, the importance of this ratio 
lies in the fact that it is related to the problem of the origin 
of the comets and that it can tell us something about the conditions 
that prevailed in the regions where these bodies were formed. 



The radical in comets affords another embarrassing example. 
I should just like to report briefly here on the work of Sauval 
(Uccle Observatory) who is studying the excitation of the X 4050 
band of the 'll - 'z electronic transition in Comet Ikeya (19631). 
Computations have been made under the simplified model assumptions 
for the (R-d and for the 'g’-d case and it has been found impossible 
to reproduce the observed intensity profile satisfactorily. The C^ 
radical having no permanent dipole moment, one would expect a 
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WAVELENGTH 


1 2 1 T ° 

Fig. 13 - High-resolution (0.15 A) scans of the C C + NH 2 4745 A 
blend in Comet Kohoutek (1973X11). The similarity of the 
comet spectrum to the laboratory spectrum is apparent. 

The two NH_ lines at the red end of the complex are 

^ 12 13 

completely beyond the shor tward-degraded C C (1,0) band 

d. 



NORMALIZED INTENSITY 


T 



4744.0 4744.5 4745.0 

WAVELENGTH IN A - 


Fig. 14 - Observed (crosses) and synthesized (solid line) 4745 A 

12 13 

blend of C C + NH 2 at high resolution. The position of 
the NH 2 lines as well as their adopted relative intensities 
are indicated; notice by comparing with the laboratory 
intensities shown in Fig. 13 that theX4744.46 line was 
appreciably weakened by- the effects of the solar absorption 
lines in this case. Also shown is the profile (dashed line) 
that would correspond to pure NH 2 emissions. 



rather high excitation temperature, similar to that of C 2 , if 
collisions are negligible. However, the observed excitation 
temperature is certainly smaller than SOO^K. In order to obtain 
theoretical spectra that bear a little resemblance to the observed 
spectra it is necessary to adopt a value for the quadrupole moment 
of which is some 10 times larger than the quadrupole moment of 
the similar molecule CO^ for instance. On the other hand, if 
collisions are important, then the exciting molecules must have a 
rather low temperature and the cross-sections for collisional 
excitation of must be much larger than for C 2 (the latter 
statement, in fact, holds not only for but also for all the 

other cometary radicals). Further study of the emissions in 
comets is certainly desirable and for this we badly need data on 
the transition probabilities for the various electronic, vibrational 
and rotational transitions. 

+ 

CH . 

I have included this molecule because it is the first ion for 

which cometary spectra have been calculated. The work here is due 

to Vreux (Li^ge) • The spectrum is very simple, being a ^ II - ' E 

transition reduced to only a few lines of the (0,0) band. Figure 15 

shows a comparison between observed and computed spectra for Comet 

Bennett (1970II). The observed spectrum is at the upper left and 

the computed spectra correspond to pure fluorescence. The agreement 

-3 

IS rather good for Fqq = 8 x 10 and p = 1.5 Debye, a value close 
to the expected value for the dipole moment of CH^ . In Figure 16 
we have again the observed spectrum and now ^-d spectra for T = 100°K (upper 
right) and T = 150'*K (lower left). These do not fit with the observed 

spectra and if we decreased the temperature below 100°K we should 
increase the ratio R(0)/R(1) and the situation would get worse. 
Whereas if we increased the temperature above 150°K, we should 
increase the ratios Q(2)/Q(l) and Q(3)/Q(2), which would also be 
going in the wrong, direction. The conclusion to be drawn from this 
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• • - • • 9 -3 

IS that the density of possible colliding molecules is << 10 cm 

4 

at a distance of about 3 x 10 km from the nucleus, where the 
spectrum was recorded. On the other hand, it may be mentioned that 
there is a puzzling situation concerning the oscillator strength 

- 3 

fQQ for which the value found here (8 x 10 ) is rather different 

from the experimental value derived from lifetime measurements (68 x 

- 3 

10 ; Smith, 1971) as well as from the upper limit obtained from 

the negative result of a search for CH^ lines in the solar spectrum, 
which gives fQQ < 1 x 10 (Grevesse and Sauval, 1971). Another 

-3 

laboratory value reported recently (11.5 x 10 ; Brzozowski et al, 

1974) is closer to the "cometary" value; it is also in better 

~3 

agreement with a theoretical estimate (6.5 x 10 ) based on 

calculations by Yoshimine et al (1973). This f-value is also of 
prime importance for the study of the interstellar medium where 
one needs to know the column densities of and of CH in 

interstellar clouds in the discussion of the formation mechanisms 
of these radicals. 

\ 

OH . 

Let us change the subject a little and consider the OH 18 cm 

radio lines in comets. I shall be very brief on this since a review 

of this question will be presented by E. Gerard shortly. The 

2 2 

fluorescence of the OH radical in the A Z - X n system has been 
analysed by Mies (1974) to obtain the populat ions of the A-sublevels 
of the ground rotational level, in which almost all OH radicals 
are, and he found that the 18 cm lines can appear either in emission 
or in absorption according to the heliocentric radial velocity of 
the comets. Then Mies has proceeded to evaluate the degree of 
polarization of the radio and uv radiations under fluorescence 
conditions. For this one must compute the population distribution 

in the 16 magnetic hyperfine states corresponding to the ground 

2 ..... 

level ^3/2* ~ 3/2. Now, provided the exciting radiation is not 

polarized, the populations of two of these states differing only 
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by the sign of the magnetic quantum number will be the same, so that 
the molecules will be aligned by solar pumping and this will produce 
some linear polarization of the radiation, the amount of polarization 
observed from the earth depending upon the phase angle. The 
maximum linear polarizations for the main radio lines (2-2 and 1-1 
transitions) are -15 % and -10 % for the absorption case observed 
in November 1973 in Comet Kohoutek (1973X11), while for the emission 
observed in January 1974 the maximum linear polarizations are - 3.5 % 
and -0.9 %. These polarizations are too small to be detected with 
present-day techniques. For the uv radiation the maximum polarization 
is 5 % and this differs by only 1.5 % from the degree of polarization 
predicted for non aligned, equilibrium populations of the hyperfine 
states. The calculations have also shown that it was a very good 
approximation to ignore the details of the distribution of the 
radicals among the hyperfine states when computing the inversion 
ratio corresponding to the two A-sublevels of the ground level. 


Let us conclude with a very brief historical note. If one were 
to describe the evolution of the study of comet spectra since the 
early days of Secchi and Huggins, he could distinguish the following 
periods : 

- the beginning were of course characterized by the endeavour to 
specify the carriers of the various emissions detected in the 
cometary spectrum; thus the now well-known radicals and atoms 
were gradually identified; this pure identification phase is in 
fact still going on, as shown by the recent discovery of H 20 ^and 
of H 2 O for instance, and further efforts should be made yet to improve 
the identification lists; 

the qualitative interpretation period started in 1911 when K. 

S chwar zs chi Id put forward the idea of fluorescence, which was 
taken again by Wurm in particular in the 1930's; this phase 
culminated in 1941 when Swings resolved the puzzle of the very 
irregular and unusual intensity profiles of the cometary molecular 
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bands by invoking the effect of the Fraunhofer lines; 
this was followed by a semi -quant i tat i ve interpretation phase 
marked by the works of McKellar and of Hunaerts in the 1940's 
and 1950's; these authors already had some success in reproducing 
the general characteristics of the cometary emissions, although 
they avoided the problem of solving the statistical equilibrium; 
the first quantitative analyses were carried out in the early 
60's soon after the first high-dispersion spectra had been secured 
by Greenstein; these studies, which are based on the solution of 
the appropriate steady-state equations, require the use of a high- 
speed computer with a large memory; 

maybe the next phase will be called the model atmosphere period 
which is still in a stage of infancy; what I have tried to show 
here is that when we enter this certainly more rewarding phase 
we are forced to be rather exacting concerning the accurate 
knowledge of all the data which are required in the calculations 
because we have to do now with secondary effects which manifest 
themselves as rather slight departures from the results obtained 
with a simplified theory; the latter point also indicates that it 
would be desirable to base these computations on observation 
material of the highest quality, possibly using photoelectric 
pho tome t ry . 
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